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We study a refrigerator model which consists of two n-level systems interacting via a pulsed 
external field. Each system couples to its own thermal bath at temperatures Th and T c , respectively 
(6 = T c /Th < 1). The refrigerator functions in two steps: thermally isolated interaction between 
the systems driven by the external field and isothermal relaxation back to equilibrium. There is 
a complementarity between the power of heat transfer from the cold bath and the efficiency: the 
latter nullifies when the former is maximized and vice versa. A reasonable compromise is achieved 
by optimizing the product of the heat-power and efficiency over the Hamiltonian of the two system. 
The efficiency is then found to be bounded from below by £ca = ,}_„ — 1 (an analogue of the 
Curzon-Ahlborn efficiency), besides being bound from above by the Carnot efficiency C,c — j^-g ~ 1- 
The lower bound is reached in the equilibrium limit 9 — ► 1. The Carnot bound is reached (for a finite 
power and a finite amount of heat transferred per cycle) for Inn > 1. If the above maximization 
is constrained by assuming homogeneous energy spectra for both systems, the efficiency is bounded 
from above by £ca and converges to it for n>l. 

PACS numbers: 05.70.Ln, 05.30.-d, 07.20.Mc, 84.60.-h 
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I. INTRODUCTION 

Thermodynamics studies principal limitations imposed 
on the performance of thermal machines, be they macro- 
scopic heat engines or refrigerators [1-3], or small devices 
in nanophysics [4] and biology [5] . Taking as an example 
a refrigerator driven by a source of work, we recall three 
basic characteristics applicable to any thermal machine: 

• Heat Q c transferred per cycle of operation from a 
cold body at temperature T c to a hot body at tem- 
perature T h (T h >T C ). 



• Power, which is the transferred heat Q c 
over the cycle duration r. 



divided 



• Efficiency (or performance coefficient) £ = Q c /W, 
which quantifies the useful output Q c over the work 
W consumed from the work-source for making the 
cycle. Note that work-consumption is obligatory, 
since the heat is transferred from cold to hot, i.e., 
against its natural gradient. 

The second law imposes the Carnot bound 

C < Cc = T c /(T h - T c ) 

on the efficiency of refrigeration [2]. Within the usual 
thermodynamics the Carnot bound (both for heat- 
engines and refrigerators) is reached only for a reversible, 
i.e., an infinitely slow process, which means it is reached 
at zero power [2, 3]. The practical value of the Carnot 
bound is frequently questioned on this ground. 

The drawback of zero power is partially cured within 
finite-time thermodynamics (FTT), which is still based 
on the quasi-equilibrium concepts [6]. For heat-engines 
FTT gives an upper bound 77 opt < ?7ca = 1 - ^/T c /T h , 



where rj op t is the efficiency at the maximal power of work- 
extraction [7] . Naturally, j]ca is smaller than the Carnot 
upper bound 1 — T c /Th for heat-engines. 

Heat engines have recently been studied within micro- 
scopic theories, where one is easily able to go beyond the 
quasi-equilibrium regime [8-15]. For certain classes of 
heat-engines the C A efficiency is a lower bound for the ef- 
ficiency at the maximal power of work [8-10]. This bound 
is reached at the quasi-equilibrium situation Th — > T c in 
agreement with the finding of FTT. The result is consis- 
tent with other studies [11, 12]. 

The interest in small-scale refrigerators is triggered 
by the importance of cooling processes for functioning 
of small devices and for displaying quantum features of 
matter [4, 15-20]. In particular, the theory of these re- 
frigerators can provide answers to several basic questions 
such as how the third law limits the performance of a 
cooling machine at low temperatures [16], and how small 
are the temperatures reachable within a finite working 
time and under a reasonable amount of resource. Natu- 
rally, the small-scale refrigerators should also operate at 
a finite power. Note that the mirror symmetry between 
heat-engines and refrigerators, which is well-known for 
the zero-power case [2], does not hold more generally [21]. 

The present situation with finite-power refrigerators is 
somewhat unclear [21-24]. Here maximizing the power 
of cooling does not lead to reasonable results, since there 
is an additional complementarity (not present for heat 
engines) [17, 18, 21, 22]: when maximizing the heat- 
transfer power one simultaneously minimizes the effi- 
ciency to zero, and vice versa. 

Here we intend to study optimal regimes of finite-power 
refrigeration via a model which can be optimized over 
almost all of its parameters. The model represents a 
junction immersed between two thermal baths at differ- 
ent temperatures and driven by an external work-source. 



This type of models is frequently studied for modelling 
heat transport; see, e.g., [4, 19, 25]. Our model is quan- 
tum, but it admits a classical interpretation, because all 
the involved density matrices will be diagonal [in the en- 
ergy representation] at initial and final moments of stud- 
ied processes 1 . 

This paper is organized as follows. The model is in- 
troduced in section II. Here we also show that the effi- 
ciency of the model is bounded by the Carnot value, and 
provide a general discussion of the refrigeration power. 
We confirm the heat-power-efficiency complementarity in 
section III and conclude that the most meaningful way 
of optimizing its functioning is to maximize the product 
of efficiency and the heat power. The optimization pro- 
cedure is reported in section IV. We discuss the quasi- 
equilibrium limit of our model in section V. There we 
show that there is a lower bound £ca = — 1 + 1/vT 
(9 = T c /Th) for the efficiency, in addition to the upper 
Carnot bound Cc = TZTr The same expression £ca was 
obtained within finite-time thermodynamics as an upper 
bound when optimizing the product of heat-power and 
efficiency or the ratio of the efficiency over the cycle time 
[21, 22]. Section VI discusses the attainability of the 
Carnot efficiency at a finite power. Entropy production 
inherent in the functioning of the model refrigerator is 
studied in section VII, while in section VIII we outline 
consequences of constraining features of the model to the 
quasi-classical domain. This constraint allows to repro- 
duce the prediction of FTT on the upper bound of Cca- 
We summarize in section IX. Some technical questions 
are relegated to Appendix. 



II. THE MODEL 

Consider two quantum systems H and C with Hamil- 
tonians Hn and He, respectively. Each system has n 
energy levels. H and C constitute the working medium 
of our refrigerator; see Fig. 1. 

Initially, H and C do not interact and are in equilib- 
rium at temperatures Th = 1/Ph > T c = 1/0 C [we set 
k B = 1]: 



P 



= e-^ Hii /tr [e~ 0hH % a = e~^ Ho /tv [e- ' H °], (1) 



where p and a are the initial Gibbsian density matrices 



1 This aspect is similar to the Ising model. This is a model for 
quantum- mechanical spin-i but it can be given a classical in- 
terpretation via an ovcrdamped particle moving in a asymmetric 
double-well potential. If the transversal components of the quan- 
tum spin are excited, this analogy breaks down. However, it still 
holds for the spin-flip process, where the transversal components 
are absent both initially and finally. In fact, the dynamics of the 
Ising model is introduced via such spin-flip processes, and this 
dynamics admits a classical interpretation. 




FIG. 1: The refrigerator model. Two systems H and C op- 
erate between two baths at temperatures T c < Th and are 
driven by an external potential V(t). W and Q c and Qh are, 
respectively, the work put into the overall system and the 
heats transferred from the cold bath and to the hot bath. 



9 of H and C, respectively. We write 



p = diag[r„,...,ri], a = diag[s„, ...,si], (2) 

Ha = diag[e„,...,ei = 0], H c = diag[/x„, ...,pi = 0], 

(3) 

where diag[a,..,6] is a diagonal matrix with entries 
(a,..., 6), and where without loss of generality we have 
nullified the lowest energy level of both H and C. Thus 
the overall initial density matrix is 



Q- m = p® a, 



(4) 



and the initial Hamiltonian Hu (gi 1 + 1 ® He ■ 

The goal of any refrigerator is to transfer heat from the 
cooler bath to the hotter one at the expense of consuming 
work from an external source. The present refrigerator 
model functions in two steps: thermally isolated work- 
consumption and isothermal relaxation; see Fig. 1. Let 
us describe these steps in detail. 

1. H and C interact with each other and with the 
external sources of work. The overall interaction is de- 
scribed via a time-dependent potential V (t, 5) in the total 
Hamiltonian 



H{t, 5) = H H <E)1 + 1(Z>H C + V(t, 6) 



(5) 



of H + C. The interaction process is thermally isolated: 
V(t, 5) is non-zero only in a short time-window < t < 
S and is so large there that the influence of all other 
couplings [e.g., couplings to the baths] can be neglected 
[pulsed regime]. The time-dependent potential V(t, 6) 
may explicitly depend on the coupling time 6. 

Thus the dynamics of H + C is unitary for < t < 6: 

n { = n(S) =UQiU\ U = Te-^o dsV{s - s \ (6) 

where f2; = 17(0) = p<g)<j is the initial state defined in (1), 
ilf is the final density matrix, U is the unitary evolution 
operator, and where T is the time-ordering operator. The 
work put into H + C reads [1, 2] 

W = Ei - Ei = tr[ {H H <8> 1 + 1 <8> H c ) (Of - fli) ], (7) 



where E{ and E[ are initial and final energies of H + C . 
2. Once the overall system H + C arrives at the final 
state Ofi n , V(t, S) is switched off, and H and C (within 
some relaxation time) return back to their initial states 
(1) under influence of the hot and cold thermal baths, 
respectively. Thus the cycle is complete and can be re- 
peated again. Because the energy is conserved during 
the relaxation, the hot bath gets an amount of heat Qh, 
while the cold bath gives up the amount of heat Q c : 

Qh = tr(ff H [tr c fif - p]), Qc = tr(H c [a - tr H Of]), (8) 

where trjj and trc are the partial traces. Eq. (1) and the 
unitarity of IA lead to 

PhQh - PcQc - S(n t ||fti) = tr[fi f lnfif - n t lnfti], (9) 

where S'(J7f||f2i) > is the relative entropy, which em- 
ployed in deriving thermodynamic bounds since [1, 35]. 

S^fif ||f2;) nullifies if and only if Of = O^ ; otherwise 
it is positive. Eq. (9) is the Clausius inequality, with 
S'(f2f||f2i) > quantifying the entropy production. This 
point will be re-addressed and confirmed in section VII. 

Eqs. (7-9) and the energy conservation Q h — W + Q c 
imply 



(fi c - h)Qc < PhW, 



(10) 



meaning that in the refrigeration regime we have Q c > 
and thus W > 0. Thus within the step 1 the work source 
transfers some energy from C to H, while in the step 2, 
C recovers this energy from the cold bath thereby cooling 
it and closing the cycle. 

Eq. (9) leads to the Carnot bound for the efficiency £ 
[we denote 6 = T c /T h < 1] 



*> W 



1 



g(Ofjja) 

(/3 C - p h )W 



< 



1 



Cc(n) 



We note from (11) that the deviation from the Carnot 
bound is controlled by the ratio of the entropy production 
S(ft f ||f2i) to the work W. 

We note in passing that all quantities introduced so far 
are meaningful also without the stage 2. Then the prob- 
lem reduces to cooling the initially equilibrium system C 
with help of the work-source and the system H. Both 
the work-source and H are clearly necessary to achieve 
cooling 2 . Q c quantifies the amount of cooling, while £ 
accounts for the relative effort of cooling. 



2 Indeed, if C and the work-source form a closed system, no cooling 
is possible due to the Thomson's formulation of the second law 
[13] (cyclic processes cannot lead to work-extraction). If H and 
C form a closed system, then W = and no cooling is possible 
due to (10). 



A. Power 

Recall that the power of refrigeration Q c /t is defined 
as the ratio of the transferred heat Q c to the cycle dura- 
tion r. For our model r is limited mainly by the duration 
of the second stage, i.e., r should be larger than the re- 
laxation time T rc i, which depends on the concrete physics 
of the system-bath coupling. 

Though some aspects of the following discussion are 
rather general, it will be useful to have in mind a con- 
crete relaxation scenario. Consider the collisional relax- 
ation scenario, where the target system interacts with 
independent bath particles via successive collisions; see 
[35, 36] and Appendix A. For our purpose the target sys- 
tem is H or C that interact with, respectively, the hot 
and cold bath. Each collision lasts a time r co i, which is 
much smaller than the characteristic time Tbtw between 
two collisions. The interaction Hamiltonian between the 
target system and a bath particle is conserved, so that 
no work is done in switching the system-bath interaction 
on and off; see [35, 36] and Appendix Al. 

The relaxation process is typically (but not always) 
exponential with the characteristic relaxation time de- 
pending on the collisional interaction; see Appendix A 2. 
This time can be much smaller than any characteristic 
time of H or C. Since the two baths act on H and C in- 
dependently, the overall relaxation process drives H + C 
to the initial state (1, 4). 

If the interaction time 5 of V(t, S) [see (6)] is also much 
smaller than Tbtw, one realizes a thermally-isolated pro- 
cess, because the overlap between the pulse and a colli- 
sion can be neglected 3 . 

To achieve a cyclic process within the exponential re- 
laxation with the relaxation time r re i, the cycle time r 
should be larger than t vc \. For each cycle the deviation 
of the post-relaxation state from the exact equilibrium 
state (1, 4) will be of order e~ T / Tro1 . Thus if the ratio 
t/t to i is simply large, but finite, one can perform roughly 
~ e T ' TTal 3> 1 number of cycles at a finite power, before 
deviations from cyclicity would accumulate and the re- 
frigerator will need resetting. 

Though, as we stressed above, the relaxation process is 



3 Analogous conclusion on the irrelevance of the system- bath inter- 
action during the pulse action is obtained when this interaction 
is always on, but its magnitude is small [weak-coupling]. Now 
the relaxation time is much larger than the internal characteristic 
time of H and C . Because the system-bath interaction is always 
on, there will be a contribution in the work (7) coming from 
the system-bath interaction Hamiltonian [37]. This contribution 
arises even when the conditions for the pulsed regime hold [37]. 
However, within the weak coupling assumption this additional 
contribution is proportional to the square of the system-bath in- 
teraction constant and can be neglected [37]. We stress that this 
additional contribution does not arise within the collisional relax- 
ation scenario, because the pulse and collisions are well-separated 
in time. 



normally exponential 4 , there are also situations within 
the collision relaxation scenario, where the system settles 
in the equilibrium state after just one inter-collision time 
Tbtw! see section VIA and Appendix A 2 for details. The 
above limitations on the number of cycles does not apply 
to this relaxation scenario. 



1. Comparing with the power of the Carnot cycle 

The above situation does differ from the power consid- 
eration of usual (reversible) thermodynamic cycles, e.g., 
the Carnot cycle [2, 3, 26, 27]. There the external fields 
driving the working medium through various stages have 
to be much slower than the relaxation to the momentary 
equilibrium. The latter means that the working medium 
is described by its equilibrium Gibbs distribution with 
time-dependent parameters. The condition of momen- 
tary equilibrium for the working medium is necessary for 
the Carnot cycle to reach the Carnot efficiency [2, 3]. 

The precise meaning of the external fields being slow 
is important here. If rp is the characteristic time of the 
fields, then the deviations from the momentary equilib- 
rium are of order Of 1 ^] [3, 26, 27]. This fact is rather 
general and does not depend on details of the system and 
of the studied process, e.g., it does not depend whether 
the process is thermally isolated or adiabatic 5 . In par- 
ticular, it is this deviation of the state from the momen- 
tary equilibrium that brings in the entropy production 

(or work dissipation) of order of 0[( ^P 1 ) ] [3, 26, 27]. 

Thus performing the reversible Carnot cycle with (ap- 
proximately) the Carnot efficiency means keeping the ra- 
tio IlsL very small. 

Now there are two basic differences between the Carnot 
cycle and our situation: 

• In our case we do not require the working medium 
to be close to its momentary equilibrium state dur- 
ing the whole process. It suffices that the medium 
gets enough time to relax to its final equilibrium. 



4 More generally, the relaxation need not be exponential, but it 
still can be such that although the difference between the system 
density-matrix at time t and the corresponding Gibbsian density 
matrix goes to zero for t — >■ oo, this difference does not turn to 
zero after any finite t. One of referees of this paper pointed out 
to us that i) the latter feature holds for a rather general class of 
relaxation processes taking place under a constant Hamiltonian; 
%%) it is rooted in the Kubo-Martin-Schwingcr (KMS) condition 
[38] for correlation function evaluated over an equilibrium state; 
see [3] for a heuristic version of this argument; Hi) the collisional 
relaxation is different in this respect, because the Hamiltonian 
is not constant. As we stress in Appendix A, a general point 
of the collisional relaxation is that no work is involved in this 
time-dependence. 

5 If a slow thermally isolated process is performed on a finite sys- 
tem, there are additional limitations in achieving the momentary 
equilibrium; sec [26, 28] for more details. These limitations are 
however not essential for the present argument. 



A small, but finite IlsL for the Carnot cycle situ- 
ation means that deviations from the momentary 
equilibrium are visible already within one cycle. 
In contrast, a small, but finite ^ sL for our situa- 
tion means that we can perform an exponentially 
large number of cycles before deviations from the 
cyclicity will be sizable. Here is a numerical ex- 
ample. Assume that lui. = I&L = 1/20. For 
the standard Carnot cycle already within one cycle 
the deviation from the momentary equilibrium will 
amount to 0.05. In our situation the same amount 
e~ 3 = 0.0498 of deviation from the cyclicity will 
come after 



= 17 _ 



2.4 x I0 7 cycles. This is a large 
number, especially taking into account that no real- 
istic machine is supposed to work indefinitely long. 
Such machines do need resetting or repairing. The 
point is that our machine can perform many cycles 
at a finite power before any resetting is necessary. 



III. COMPLEMENTARITY BETWEEN THE 
TRANSFERRED HEAT AND EFFICIENCY 

We now proceed to optimizing the functioning of the 
refrigerator over the three sets of available parameters: 
the energy spacings {£fc}5J =2 , {Mfc}fc = 2' ana - the un it ar y 
operators U. It should be evident from (5, I) that op- 
timizing over these parameters is equivalent to optimiz- 
ing over the full time-dependent Hamiltonian H(t, S) of 
H + C. We stress in this context that no limitations on 
the magnitude of V(t, S) are imposed. This means that 
the unitary operator can in principle be generated in an 
arbitrary short coupling time S. 

We start by maximizing the transferred heat Q c = 
tr(i?c[<7 — trn^f]), which is the main characteristics of 
the refrigerator. Since tr[ffc(r] depends only on {pk}k = 2> 
we choose {£/c}JJ =2 an< ^ ^M so that the final energy 
tr[i?c^f] attains its minimal value zero. Then we maxi- 
mize tracer] over {pk}k=2- Note from (2) 



1® -He = diag[/ii ,... 
fl[ = p<E)(J = diag[siri, . . 



, [Ai, . . . , /i n > • ■ 

) SlTfi i • • • i °n'\ i 






It is clear that tr[i?c^f] = ti:[HcUQ,\U^] goes to zero 
when, e.g., r 2 = . . . = r n ->• (e = e 2 = • • • = e n ->• oo), 
while U amounts to the SWAP operation Up ® aW = 
a ® p. It is checked by a direct inspection that the maxi- 
mization of the initial energy ti[Hc<j] over {^fe}fe =2 pro- 
duces the same structure of n — 1 times degenerate upper 
energy levels jj, = /X2 = • • ■ = jJ> n ■ Denoting 



V = S 2 = .. = i 

we obtain for Q 
Q c = T c In 



-/3c/* 



r-2 



-0 h s 



(v — u)(n — 1) 



l + (n- l)v][l + (n - l)u] 



(12) 



(13) 



where according to the above discussion, Q c is maximized 
for u — > 0, and where v is to be found from maximizing 
Qc|«->o m (13) over v, i.e., v is determined via 



l + (»- l)v + lnw = 0. 



(14) 



For the efficiency we get for the present situation (H 
and C have n—1 times degenerate upper levels, while hi 
amounts to the SWAP operation): 



Qc 
W 



/' 



OMW 



e — fi 



Hi 



01n[ 



(15) 



The maximization of Q c led us to u — > 0, which then 
means that £ in (15) goes to zero. 

Thus C can be cooled down to its ground state 
(tr[i?cOf] — >• 0), but at a vanishing efficiency, i.e., at 
expense of an infinite work. To make this result consis- 
tent with the classic message of the third law [33], we 
should slightly adjust the latter: one cannot reach the 
zero temperature [of an initially equilibrium system] in a 
finite time and with finite resources [infinite work is not 
a finite resource]. At any rate, one should note that the 
classic formulation of the third law motivates its opera- 
tional statement using exclusively equilibrium concepts. 
Modern perspectives on the third law are discussed in 
[16, 18, 20, 29, 30]. 

Note that the efficiency £ in (15) reaches its maximal 
Carnot value 0/(1 — 9) for 



u = v, 



(16) 



which nullifies the transferred heat Q c ; see (13). 

Now we have to show that Q c tends to zero upon max- 
imizing £ over all free parameters {£fc}£ =2 , {pk}k=2 an d 
U. Denoting {|«H)}fc = i and {|«c)}JJ = i for the eigenvec- 
tors of Hh and He, respectively, we note from (7, 8) that 
W and Q c feel U only via the matrix 



a 



ij | kl 



KinjcMknlc) 



This matrix is double-stochastic [34] 



J2 H c ^\ki-J2 kl a ' 



ij | kl 



= 1. 



(17) 



(18) 



Conversely, for any double-stochastic matrix Cy | ki there 
is some unitary matrix U with matrix elements U^j i ^ , 



so that C, 



ij I kl 



= \u ij]kl \ 2 



[34]. Thus, when maximizing 
various functions of W and Q c over the unitary hi, we 
can directly maximize over the (n 2 — l) 2 independent 
elements of n 2 x n 2 double stochastic matrix CV,- 1 f-i . 

We did not find an analytic way of carrying out the 
complete maximization of £ over all free parameters. 
Thus we had to rely on numerical recipes of Mathemat- 
ica 7, which for n = 1, . . . , 5 confirmed that Q c nullifies 
whenever £ reaches (along any path) its maximal Carnot 
value. We believe this holds for an arbitrary n, though 
we lack any rigorous proof of this assertion. 




FIG. 2: Solid line: efficiency ( of the optimized refrigerator 
versus the temperature ratio 8 = T c /Th for n = 3; see (15). 
In the scale of this figure f(n = 2) and f (n = 3) are almost 
indistinguishable. Dashed line: the lower bound J= — 1. 



IV. MAXIMIZING THE PRODUCT OF THE 
TRANSFERRED HEAT AND EFFICIENCY. 

We saw above that neither Q c nor £ are good target 
quantities for determining an optimal regime of refriger- 
ation. But 



X : 



cC, 



(19) 



is such a target quantity, as will be seen shortly. This is 
the most natural choice for our setup. This choice was 
also employed in [21]. Refs. [17, 18, 24] report on different 
approaches to defining refrigeration regimes. 

The numerical maximization of \ — QQc over {£fc}5J =2 , 
{fJ-k}k=2 an d U nas been carried out for n — 1,...,5 
along the lines discussed around (17, 18). It produced 
the same structure: both H and C have n—1 times 
degenerate upper levels, see (12), and the optimal U again 
corresponds to SWAP operation 6 : 



fii = p a, ilf = 



(20) 



Recalling the expression (5) for the total Hamiltonian 
H(t, 6), we can state this result as follows: there exist a 
coupling potential V(t, 8) that for a given coupling time 
8 generates the unitary SWAP operation following to (6). 
This operation does not explicitly depend on <5, because 
V(t, 8) depends on 5; see also our discussion in the begin- 
ning of section III. Note that both the initial and final 
states in (20) are diagonal in the energy representation. 



6 Let us recall how the SWAP is defined via a pure-state base. Let 
{|^')i}fc=i be an ortho-normal base in the Hilbert where p lives. 
Let also {|A;)2}J!_ 1 be an orthonormal base in the Hilbert space 
where a lives. Any unitary operator acting on the composite 
Hilbert space can be defined with respect to the orthonormal 
base |fc)i ® IO21 where fe, I = 1, . . . ,n. Let us now define for all 
pairs k and /: !7swAp|fc)l ® IO2 = \l)l ® |fe)2. 



Evidently, the intermediate state fl(t) for < t < 5 is 
not diagonal in this representation. 

The efficiency C, and the transferred heat Q c are given 
by respectively, (13) and (15), where instead of u and 
v we should substitute u and v, respectively. The latter 
two quantities are obtained from maximizing \ — QQ C , 



X(u,v) 



T c 0{n- l)(v-u)ln 



2 1 



[In I _ 0]n I][l + (n - !)«][! + (n - l)v] 



,(21) 



where u and v are found from maximizing x( u j v ) Y ' ia - 
duX = d v x = 0. Note that u and t> depend on 9 = T c /Th- 
The efficiency C, and the transferred heat Q c are given, 
respectively, by (15) and (13) with u — > u and v — > v. 

Though we have numerically checked these results for 
n < 5 only, we again trust that they hold for an arbitrary 
n (one can, of course, always consider the above structure 
of energy spacings and hi as a useful ansatz). 

SWAP is one of the basic gates of quantum informa- 
tion processing [31]; see [14] for an interesting discussion 
on the computational power of thermodynamic processes. 
SWAP is sometimes realized as a composition of more el- 
ementary unitary operations, but its direct realizations 
in realistic systems also attracted attention; see, e.g., [32] 
for a direct implementation of SWAP in quantum optics. 
Note that for implementing the SWAP as in (20) the ex- 
ternal agent need not have any information on the actual 
density matrices p and a. 



A. Effective temperatures 

Since the state f2f of H + C after the action of V(t) 
is a p, and because in the optimal regime the upper 
level for both H and C is n — 1 times degenerate, one 
can introduce non-equilibrium temperatures T' h and T' c 
for respectively H and C via [note (1)] 



a = e 



-^ H "/tv[e-^ H "], p^e-^ Hc /tT[e-^ Hc ],(22) 



where we recall that a (p) is the state of H (C) after 
applying the pulse. Using (12) we deduce 

lni , lni 

V u 

where v = e~^ c/1 and u = e~^ h£ ; see (12). This implies 

T c T h = T' c T' h . (24) 

As expected, the refrigeration condition v > u, see (13, 
21), is equivalent to 



T' c <T c <T h < Ti 



(25) 



i.e., after the pulse the cold system gets colder, while 
the hot system gets hotter. Note that the existence of 
temperatures T' c and T' h was not imposed, they emerged 



out of optimization. In terms of these temperatures the 
efficiency (15) is conveniently written as 



c 



T r 



V, 



T' h - T c T h - T' c 



(26) 



We eventually focus on two important limits: quasi- 
equilibrium 6 — ¥ 1, and the regime lnn> 1. 



V. QUASI-EQUILIBRIUM REGIME 9 -»• 1: A 
LOWER BOUND FOR THE EFFICIENCY 



In this regime the temperatures T^ and T c are nearly 
equal to each other: = T c /Th — *• 1. 

First we note that sharply at $ = 1, \ reads 

x(o)|fl=i = T c 6{n - 1) [1 + (n - l)a]" 2 a In 2 a, (27) 

where 

u ~ v = a, 

and where a is given by d a x{a)\o=i = 0: 

[(n-l)o-l] lna = 2[(n-l)o+l]. (28) 

We now work out the optimal u and v for 9 — > 1. It can 
be seen from (21) that the proper expansion parameter 
for 9 — *• 1 is x = Vl — 0. We write 

u = a + y^ a k x k , v = a + y^ (a k + b k ^ 1 )x k .(29) 

-*- — 4 k— 1 -*- — J k— 1 

We substitute (29) into d u x = and d v x = and expand 
them over x. Both expansions start from terms of order 
O(x ). Now a k and b k are determined by equating to 
zero the 0(x k ) terms in d u x — and d v \ = 0. Thus the 
O(x ) terms together with (28) define 6q : 



bo = a In — , 



(30) 



which should be non- negative due to v > u. The 0(x l ) 
terms together with (28) and (30) define a\ and b±: 

a , 1 , a[24 + In 2 a] , 1 
oi = --In-, 6 1 = — L-I 1 In-, 31 

2 a 48 a 

and so on. Eqs. (29, 30, 31) imply for the efficiency at 
6 -> 1 (x = Vl~9) 

1 In 2 a [48 + In 2 a] In 2 a , 9n , s 

C — - 1 + ^8~- 1536 X + °^- ^ 

Note that the expansion (32) does not apply for n — > oo, 
since in this limit a(n — 1) ~ 1 — ln i*_^ ', see (28). Thus, 

in the limit 9 — > 1, Q c scales as ex \f\ — 9, 
aT c [\na] 2 (n-l)VT~9 



;i + (n-l)o] 



(33) 



while the consumed work is smaller and scales as 1 — 9. 
Eq. (32) suggests that the maximization of \ imposes 
a lower bound on the efficiency: 



C>Cca 



f 



VT 



1. 



(34) 



This is numerically checked to be the case for all < 9 < 
1 and all n; see also Fig. 2. 

The expression of £ca was already obtained within 
finite-time thermodynamics — but as an upper bound on 
the efficiency — and argued to be an analogue of the 
Curzon-Ahlborn efficiency for refrigerators [21, 22]. Sec- 
tion VIII explains that also within the present micro- 
scopic approach £ca can be an upper bound for £ pro- 
vided that x i s maximized under certain constraints. 

Recalling (15), our discussion after (23-25) and (26), 
we can interpret the lower bound for the efficiency as a 
lower bound on the intermediate temperature T' of C: 



1-VT=- 



T. 
< — < 1, 



(35) 



i.e., the lowest temperature T' c cannot be too low un- 
der optimal x- Compare this with the fact that under 
vanishing efficiency (that is for very large amount of the 
consumed work), T' c can be arbitrary low; see our dis- 
cussion after (15). Thus a well-defined lowest (per cycle) 
temperature emerged once we restricted the resource of 
cooling (the consumed work). 



VI. THE MANY-LEVEL REGIME: REACHING 
THE CARNOT LIMIT AT A FINITE POWER. 



Now we turn to studying the regime 
ln(n - 1) > 1. 
First of all, let us introduce two new variables 

_ ln[n - 1] 



denote 



u(n — 1) ln[n — 1], £ : 



p=ln[ra- 1], 



v (n — 1) ' 



and rewrite x m (21) as 



x(i-fl 

9T cP 



(1 



S)(l 



±ln[£]) s 



(i + f)(i + |)(i + ir%M&) 



(36) 



(37) 



(38) 



(39) 



The expression in RHS of (39) is now to be optimized 
over p and £. We note that if these parameters stay 
finite in the limit p = ln[n — 1] 3> 1, the value of x is read 
off directly: g T ~ ' = 1. The finitness of p and £ in the 
limit p = ln[n — 1] 3> 1 is confirmed by expanding the 
RHS of (39) over the small parameter -, collecting terms 



0.105 - 



0.100 - 



0.090 




FIG. 3: Convergence of the efficiency £ (normal line) to the 
Carnot value i]c — 1/9 (dashed line) as a function of i^- 2 , 
where p = ln[n — 1]; see (36-39). 

oc O(-), differentiating them over p and £, and equating 
the resulting expressions to zero. This produces: 



P = 



1- 



O 



i = 



2- 



o 



Substituting these into (15) and (21) we get 



Qc 

T, 



1 -i 
3 



26» ln[p] 
{l-9) 2 ~p~ 



l-i 



In 



2-i 



O 



(40) 

(41) 
(42) 



Note from (36, 38) that the dominant factor in the ef- 
ficiency C is the Carnot value tzb, while the subleading 
term is naturally negative; see (41). Likewise, the dom- 
inant factor in Ur- is p = ln[n — 1], while the subleading 
term is 0(1). We also see that the limit ln(n — 1) 3> 1 
does not commute with the equilibrium limit 9 — > 1, since 
the corrections in (41, 42) diverge for 6 — > 1. 

Thus in this regime p = ha.(n — 1)3>1 the efficiency 
converges to the Carnot value; see Fig. 3. 

Recalling (12) we see from (37, 40) that in the regime 
(36) the total occupation of the higher levels of H is 
small, so that H is predominantly in its ground state 
before applying the work-consuming external pulse. In 
contrast, C is more probably in one of its excited states. 
These facts are expected, because C has to give up some 
energy, while H has to accept it. 

We note that this regime resembles in several aspects 
the macroscopic regime N 3> 1 of a ./V-particle system. 
Recall that for N ^> 1 (weakly coupled) particles the 
number of energy levels scales as e N , while energy scales 
as N. Now for the above situation (42) the transferred 
heat Q c is (in the leading order) a product of the colder 
temperature T c and the "number of particles" ln(n — 1). 

The effective temperatures T' h and T' c [see (23-25)] in 



this limit are close to their initial values: 



rpl 


= 1 + - In 
P 


"(i-0)V 

2-6 


+ o 


' 1 " 

p 2 


(43) 


V 

■*■ C 


= 1 - - In 
P 


"(1-0)V 
2-6 


+ o 


' 1 ' 

p 2 


(44) 



where we employed (37, 40). Though during the refrig- 
eration process the systems C and H are able to process 
large amounts of work and heat (oc ln[n — 1]), their tem- 
peratures are not perturbed strongly. 



A. Finiteness of power 

It is important to note that the asymptotic attainabil- 
ity (41) of the Carnot bound for £ is related to a finite 
transferred heat Q c = T c ln[n — 1], but it also can be 
related to a finite power -^, where in our model the cy- 
cle time r basically coincides with the relaxation time; 
recall our discussion in section II A. This appears to 
be unexpected, because within the standard thermody- 
namic analysis the Carnot efficiency is reached by the 
Carnot cycle at a vanishing power [2] ; see section II A 1 
for a precise meaning of this statement. In any refrig- 
erator model known to us — see, e.g., [19] — approaching 
the Carnot limit means nullifying the power. See also in 
this context our discussion around equation (16); various 
reasons preventing the approach to the Carnot efficiency 
for thermal machines (even for small machines working at 
zero power) are analysed in [26] . Now we supplement our 
discussion in section II A with more specific arguments. 

We already stressed in section II A that within the sec- 
ond stage of the refrigerator functioning, where both H 
and C relax to equilibrium under influence of the corre- 
sponding thermal baths, the relaxation mechanism can 
be associated with the collisional system-bath interac- 
tion; see Appendix A for a detailed discussion of this 
mechanism. Here there are three characteristic times: 
the single collision duration time r co i is much smaller 
than the inter-collision time Tbtw, while the relaxation 
of the system to its equilibrium state is governed by the 
time r re i. The assumed condition t co \ <C Tbtw allows im- 
plementing the thermally isolated work-consuming pulse, 
because if the pulse time is also much smaller than Tbtw, 
the pulse does not overlap with collisions. 

In Appendix A 2 we study the relaxation time of the 
system with n — 1 fold degenerate upper energy levels 
and non-degenerate lowest energy level. We also account 
for the limit ln[n— 1] 3> 1, where the Carnot efficiency is 
reached; see (41). It is shown that for such a system the 
relaxation time r ro i can — depending on the details of the 
thermal bath and its interaction with the system — range 
from few Tbtw's to a very long times oc riThtw The former 
relaxation time means a finite power, while the latter 
time implies vanishing power oc ^-^ for lnn> 1. 

These two extreme cases are easy to describe without 
addressing the formalism of Appendices Al and A 2. 



1. For simplicity let us focus on the relaxation of the 
system H that after the work-extracting pulse (20) is 
left in the state a, and is now subjected to a stream 
of the bath particles (the situation with C is very simi- 
lar). Recall that each bath particle before colliding with 
H is in the Gibbsian equilibrium state at the tempera- 
ture Th- Now assume that each bath particle also has 
n — 1 fold degenerate upper level, and one lowest energy 
level. Also, the non-zero energy spacing for the bath 
particles is equal to that of H. Then the relaxation of 
H is achieved just after one collision provided that the 
system-bath interaction [during this collision] amounts to 
a SWAP operation. Note that the characteristic time of 
this relaxation is Tbtw, and that this is a non-exponential 
scenario of relaxation, because the system exactly settles 
into its equilibrium Gibbsian state after the first collision. 

No work is done during collisional relaxation; see Ap- 
pendix A. Indeed, under above assumptions on the en- 
ergy levels of the bath particles, the SWAP operation 
commutes with the free Hamiltonian H\ 1 + 1 Hn 
(where H\ is the Hamiltonian of the given bath particle) , 
which implies that the final energy of H plus the bath 
particle is equal to its initial value. Since each separate 
collision is a thermally isolated process, this means that 
no work is done; see (7). 

2. If each collision is very weak and almost does not 
exchange heat with the system H, the relaxation time 
becomes very long. Intermediate cases are discussed in 
Appendix A 2. These intermediate cases are relevant, 
since the power of refrigeration is finite even for long 
relaxation times oc In n. Indeed, we recall from (42) that 
Q c = T c ln[n-l] + 0(l). 



B. More realistic spectra still allowing to reach the 
Carnot bound 



One can ask whether the convergence (41) to the 
Carnot bound is a unique feature of the spectra (12) in 
the limit ln[n — 1] 3> 1, or whether there are other sit- 
uations that still allow £ — > (c- The answer is positive 
as we now intend to show. For the energy spectra (3) we 
postulate [k — 1, . . . , n — 1] 

e k +i = e + (k - 1)5, n k+1 = fj, + (k - 1)6, (45) 

where <5 > is a parameter. Next, we assume that the 
following six conditions hold 

(n - 1)135 > 1, /3 C <5 < 1, (3 h 5 < 1 (46) 

p = ln[T c /<J]»l, (47) 

u = e- Ph€ ~ 5/p, v = e~ 0c ^ ~p5. (48) 

Under conditions (46, 47, 48) — and assuming the SWAP 
operation for the pulse — we show below that the results 
analogous to (41, 42) hold, 



c = 



o 



1" 


Qc 


( 


r 11 


p. 


' T~ c 


= p(i + o 


p. 



(49) 



where the role of a large parameter p = ln[n — 1] in (41, 
42) is now played by p 7 . Note that (46) and (47) still 
imply that ln[n — 1] 3> 1. 

The spectra (45) under condition (46, 47, 48) cor- 
respond to a quasi-continuous part separated from the 
ground state by a gap. This type of spectrum avoids 
the strong degeneracy of (12), and is met in conventional 
superconductors below the transition temperature [39] . 

To derive (49) via (45-48) we note the following for- 
mulas for the partition sums Z^ = tr [e^ l3hHli ] and 
Z c = tr [e-P<= Hc ] in (1), the heat Q c and work W: 



z h -- 


= i + 7 - 


S 


Qc- 


= T c ln 


1 



Z r = 1 



S ' 



W = T h 
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VII. 



ENTROPY PRODUCTION 



Entropy production is an important characteristics of 
thermal machines, because it quantifies the irreversibility 
of their functioning [1, 2]. For our refrigerator model, no 
entropy is produced during the first stage, which is ther- 
mally isolated from the baths. However, a finite amount 
of entropy is produced during the second, relaxational 
stage. The overall entropy production reads 

&i = b [S ifi n 1 1 S Jj n J , 

and controls the deviation of efficiency from its maximal 
Carnot value; see (11). In the optimal conditions (20, 
12), we get 



Si 



S[p\\a) + S[a\\p} 

\n[v/u] (v — u){n 



1) 



l + {n- l)v][l + (n- 1)- 



(50) 
(51) 



where 5[cr||p] and S[/o||ct] are the entropies produced in, 
respectively, cold and hot bath. Indeed, consider the 
system C that after the external field action is left in 
the state with density matrix p [see (20)], and now under 
influence of the thermal bath should return to its initial 
state a ex e~P cHc . Now 



T c S[p\\a] = tr[H c p] - T c S[p] + T c \n tr< 



cH c 



(52) 



is the difference between the non-equilibrium free energy 
tr[Hcp] — T c S[p] of C in the state p [and in contact to 
a thermal bath at temperature T c ] and the equilibrium 
free energy — T c ln tr e^^ cHc . Simultaneously, T c S'[/9||cr] 



7 In the definition (47) of p one can as well employ 7), instead 
of T c . This will not lead to serious changes, because we always 
assume that 8 = T c /Th is fixed. 



in (52) is the maximal work that can be extracted from 
the system C (in state a) in contact with the T c -bath [8] . 
During relaxation this potential work is let to relax into 
the T c -bath increasing its entropy by *S"[/o| | ex] . Likewise, 
5[cr||p] is the entropy production during the relaxation of 
the system H in contact with the T/j-bath. 

Now in the regime ln[n — 1] 3> 1, Si amounts to 
ln[w/u] ~ ln(ln[n— 1]), see (51), while the consumed work 
W and the transferred heat Q c scale as ln[n— 1]. In other 
words, the entropy production Si is much smaller than 
both W and Q c . This explains why for a large ln[n — 1] 
the Carnot efficiency is reached; see (11, 41). 

In the equilibrium limit 6 — >• 1, Si reads 



S,= 



[lna] 2 (n-l)(l-<9) 
[l + (n-l)a] 2 : 



(53) 



where a is given by (28), and where in deriving (53) we 
employed (51) and asymptotic expansions presented after 
(28). Note that now Si is smaller than Q c ex y/1 — 9, but 
has the same scale 1 — 9 as the consumed work W; see our 
discussion after (32). Thus Si cannot be neglected, and 
this explains why the Carnot efficiency is not reached in 
the equilibrium limit 6—^1; sec (11). 



VIII. CLASSICAL LIMIT. 

We saw above that the optimization of the target quan- 
tity x — QcC produced an inhomogeneous type of spec- 
trum, where a batch of (quasi)degenerate energy levels is 
separated from the ground state by a gap. It is meaning- 
ful to carry out the optimization of \ imposing a certain 
homogeneity in the spectra of H and C. The simplest 
situation of this type is the equidistant spectra 



(n-l)e, p n = (n-l)p, 



(54) 
0, p, -> 



for H and C; recall (3). For n — >• oo and e — )• 
these spectra correspond to the classical limit. 

Thus, now we maximize \ = QcC imposing conditions 
(54) . We found numerically that the optimal U again cor- 
responds to SWAP operation; see (20). For \ — %(u, v) 
we get 



T r 9 In 2 i 



X 



lni 



»lni 



n(v r 



(l-w)(l-u) (1 -«")(! -u n ) 



where v = e~P" e and u = e _/3h ^ are found from maxi- 
mizing x- The efficiency C is still given by (15). 
In the limit n>l we get from maximizing \: 



1. 



1 and 



n(v T: 



(l-u n )(l -u n ) 



0, (55) 



implying that \ an d C depend on one parameter <p = j^ : 
T c 0(4>-1) , 9 



10 



The optimal value of this parameter is 
This leads to 



T r > 



X = 



(i + VT~8) 2 ' 



C = CcA = 



1 



vT^fl 



1 + vr 



- 1. (56) 



Thus for a large number of equidistant energy levels the 
maximization of \ leads to homogeneity (e — > 0, /2 — > 0), 
which is an indication of the classical limit. The efficiency 
C in this constrained optimal situation is equal to Cca- 

The above results refer to optimizing \ in the limit 
n ^> 1. However, we confirmed numerically that the 
above values (56) for \ and £ — obtained in the limit n 3> 
1 — are upper bounds for \ and C at a finite n. 

Our conclusion is that the efficiency Cca — which is a 
lower bound for the efficiency during the unconstrained 
optimization of \ — appears to be an upper bound for 
refrigerators that operate under equidistant (classical) 
spectra. The upper bound is reached in the limit n^$> 1. 

These facts clarify to some extent why Cca was ob- 
tained as an upper bound for the efficiency within the 
finite-time thermodynamics (FTT) [21, 22]. Apparently, 
the quasi-equilibrium assumptions of FTT implied con- 
straints which are equivalent to imposing the homoge- 
neous spectra in our approach. 



IX. SUMMARY 

We have studied a model of a refrigerator aiming to 
understand its optimal performance at a finite cooling 
power; see Fig. 1. The structure of the model is such that 
it can be optimized over almost all its parameters; addi- 
tional constraints can and have been considered, though. 
We have confirmed the complementarity between opti- 
mizing the heat Q c transferred from the cold bath T c and 
efficiency £: maximizing one nullifies the other. Similar 
effect for different models of quantum refrigerators are 
reported in [17-19]. 

To get a balance between Q c and £ we have thus chosen 
to optimize their product C,Q C . This leads to a lower 
bound Cca = ,}_ B — 1 (6* = |f ) for the efficiency in 



n^e 



addition to the upper Carnot bound £c = jzrg ~~ 1- The 
fact of £ > Cca implies that there is the lowest finite 
temperature reachable within one cycle of refrigeration; 
see (35). 

The lower bound Cca is reached in the equilibrium 
limit T c —> Th- Constraining both systems to have ho- 
mogeneous (classical) spectra, Cca is reached as an up- 
per bound. This is just like within finite-time thermo- 
dynamics (FTT), when maximizing the product of the 
cooling-power and efficiency [21], or the ratio of the effi- 
ciency and the cycle time [22] . In this sense Cca seems to 
be universal. It may play the same role as the Curzon- 
Ahlborn efficiency for heat engines rycA, which, again, is 
an upper bound within FTT [7], but appears as a lower 
bound for the engine models studied in [8]. For other 
opinions on the Curzon-Ahlborn efficiency for refrigera- 
tors see [23, 24]. 

The Carnot upper bound is asymptotically reached in 
the many-level limit of the model. We saw that this 
asymptotic convergence is related to a finite heat trans- 
ferred per cycle, and we argued that it can also be related 
to a finite power if the relaxation scenario of the model 
refrigerator is chosen properly: provided that the cycle 
time is larger than the relaxation time, one can perform 
exponentially large number of refrigeration cycles before 
inevitable deviations from cyclicity — that in any case are 
there due to a finite cycle time — will accumulate. To our 
knowledge such an effect has never been seen so far for 
refrigerator models. 

For the optimal refrigerator the transferred heat Q c 
behaves as Q c ex T c (in particular, for T c —J- 0); see 
(13, 21, 41). This is in agreement with the optimal low- 
temperature behaviour of Q c from the viewpoint of the 
third law [16, 20]. 
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Appendix A: Collisional relaxation. 
1. General consideration 

The purpose of this discussion is to outline the general 
structure of a collisional relaxation process. Our presen- 
tation follows to [35, 36]. 

The thermal bath is modeled as a collection of N ^> 1 
independent equilibrium systems (particles) with initial 
density matrices w, = g- cxp[— f3Hi] and Hamiltonians 
Hi, where i = 1, .., N, and where 1//3 — T is the bath 
temperature. This formalizes the intuitive notion of the 
bath as a collection of many weakly-interacting particles. 

The target system H starts in [an arbitrary] initial 
state ph and has Hamiltonian Hh. The collisional re- 
laxation is realized when the particles of the bath se- 
quentially interact [collide] with H. Multiple collisions 
(between the target system and simultaneously two or 



more bath particles) are neglected. 

Consider the first collision. The initial state of H and 
the first bath particle is Oi + h = p<E)0Ji- The interaction 
between them is realized via a unitary operator V, so 
that the final state after the first collision is f^ , H = 
V£1i+hV^. This unitary operator is generated by the full 
Hamiltonian H\ + n: 



-Hi+h — H\ + Hh + Hi h, 



(Al) 



where H\ . h is the interaction Hamiltonian. Define sep- 
arate final states: 



trifii. 



Hi 



tr H ni. 



H: 



(A2) 



where tri and trn are the partial over the first parti- 
cle and H, respectively. Recall the definition (9) of the 
relative entropy. The unitarity of V implies 



S[Cl' u 



H 



Wi] = tr[fii +H lnOi +H ] 

- tr[fti +H in(//®wi)]. 



(A3) 



Employing u>\ = -%- exp[— (3Hi] and S[ ^' 1+H 1 1 p' u>i ] > 
in (A3) we get 



TAS H + AE/i > 0, 



(A4) 



where AS*h = tr [— //In// + plnp] and AUi — 
tr(Hi [gl>i — Wi]) are, respectively, the change of the en- 
tropy of R and the average energy of the first particle. 

We now require that the interaction V conserves the 
average energy: 



AE/i = ■ 
Using this in (A4) one has 



-AU, 



AU H - TAS H < 0. 



(A5) 



(A6) 



Since we did not use any special feature of the initial 
state of H, (A6) holds for subsequent collisions of H with 
the bath particles. Thus Uh — TSu decays in time, and 
it should attain its minimum. It is well-known [35, 36] 
that this minimum is reached for the Gibbs matrix p oc 
e -0ff H . collisions can drive H to equilibrium starting 
from an arbitrary state [35, 36]. 

Condition (A5) expresses the average energy conser- 
vation. It is natural to use a more stringent condition 
according to which the sum of energies of H and the 
bath particle 1 is conserved in time [36]: 



[#i+#h,#i,h]=0. 



(A7) 



This condition makes the dynamics autonomous, since 
for any initial state of H + 1 the switching the inter- 
action on and off does not cost energy and (A5) holds 
automatically. 

For condition (A7) to be non-trivial, the operator H\ + 
Hh should have a degenerate spectrum. Otherwise due 
to [Hi + Hn,Hi] = and (A7), H H and Hi will be 
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constants of motion, which means that no transfer of 
energy and thus no relaxation is possible. 

Here are two crucial points of the collisional relaxation. 

1. If the target system starts in the equilibrium state, 
this state does not change in time under subsequent 
collisions. This analogue of the zero law of thermody- 
namics is especially obvious from condition (A7), but 
it also holds simply from the conservation of the aver- 
age energy (A5). Indeed, if pu is the equilibrium state, 
Ph oc e~P HiI , and also condition (A5) holds, the relative 
entropy S[n' 1+Ii \\ pu Wi] is equal to zero, which can 
happen only for f2' 1+H = pu wi . 

2. No work is done for switching collisions on and off. 
This is clearly seen from (A7) , which states that the free 
Hamiltonian is a constant of motion. 



2. The relaxation time for a pertinent example 

Now we study the relaxation time for an n-level sys- 
tem H under collisional dynamics. We assume that n — 1 
levels of H coincide and have energy e > 0. The lowest 
energy level is not degenerate and has energy zero. Im- 
portantly, we assume that condition (A7) holds meaning 
that the relaxation proceeds autonomously, i.e., without 
additional energy [work] costs on switching the interac- 
tion on and off. The initial (before colliding with the 
first bath particle) density matrix of H is assumed to be 
Gibbsian at temperature Tq — 1/(3q: 



-l3oH r 



ZhWo) 



rPo 



1-r 
n- 1 



1 



l + (n- l)e-^ e ' 



where Po = |0)(0| and P e is the projectors on the n— 1- 
dimcnsional eigen-space of p with eigenvalue e. 

To satisfy the degeneracy of the interaction Hamilto- 
nian [see our discussion after (A7)] we assume that the 
first bath particle has (among others energies) energy lev- 
els E and E + e. The degeneracies of these levels are rvj 

and n E+e , respectively. The equilibrium density matrix 
of the bath particle 1 is written as 



Pl=P 1+ ^P™ 



■4 P ^- (A8) 



" ' * rv " 



(A9) 



E 



where rU 1 and r E+e are the Boltzmann weights for the 



energy levels E and E + e, respectively, the summation in 
(A9) is taken over all energy levels of the bath particle, 

P E and PgJ.^; are the projectors on the corresponding 
sub-spaces, 



trpW= n W, 



trP 



[i] 



i [1] 



(A10) 



and where J>\ in (A8) is the remainder of p\. 

It is assumed that the unitary operator V responsible 
for the interaction operates within the sub-space with 



the projector P £ P^ + Po Pg +e (this sub-space has 
energy E + e), i.e., 

[V,P E ®PW+P ®P™ + J=0. (All) 

Then the post-collision density matrix p' of H reads 
p' = triVp0piV f = p 
111 L ~' A [nW Po-t ri VPo0P™ Vt] 



[i] 



- r TV,. - r 



E+e ' E n _l I V-E+e 1 



Clearly, p! commutes with Hn. For simplicity, we choose 
V such that the degeneracy of p is not resolved, i.e., in 
the state p', the occupations of the higher energy levels 
of H are equal. This means we need to keep track of the 
lowest energy-level occupation (0|p'|0) = r' only: 



<• — r — —A \r — 



1 



cqj i ' cq 



1 + (n - l)e-0' 



-, (A12) 



A 



r eq (n- 1) 



i ll] 

i E+e 



(tr 1 VP o 0P^ 1 +£ V t )|O) 



(A13) 



where r cq is the equilibrium value of r. A can be maxi- 
mized over the unitary V [under condition (All)] 



4 — 

n max — 



,w 



W(n-l) 



' cq 



(n-l) 



(A14) 



Using (A9) one can show that A < A max < 1: after 
first collision H gets closer to its equilibrium state; see 
(A12). This equation obviously generalizes to subsequent 
collisions [we revert from (A15) to (A12) for m = 1]: 



[m-l] 



= (1 - A) m [r 



' cqj 



(A15) 



It is seen that (A15) predicts exponential (with respect to 
the number of collisions) relaxation towards the equilib- 
rium value r cq of r. The approach to equilibrium is gov- 
erned by the factor (1 — ^4)" meaning that when |^4| <§C 1 
the effective number of collisions after which the equilib- 
rium is established equals to — l/[ln(l — A)]. 

Now the shortest relaxation corresponds to just one 
collision and it is reached for A = 1, e.g., r E = r cq and 

n E+e = n — 1 in (A14). Then the corresponding unitary 
operator V is the SWAP operation. The relaxation time 
in this case amounts to one inter-collision time. 

It should be clear that there is no upper limit on the 
relaxation time. The latter can be arbitrary large, e.g., 
due to V converging to 1 in (A13). Various intermediate 
cases can be studied with help of (A14). In particular, it 
is not difficult to identify regimes, where the relaxation 
time scales as ex Inn. 



